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Abstract. We investigate the process of circularly polarized high harmonic 
generation in molecules using a bicircular laser field. In this context, we show that 
molecules offer a very robust framework for the production of circularly polarized 
harmonics, provided their symmetry is compatible with that of the laser field. Using 
a discrete time-dependent symmetry analysis, we show how all the features (harmonic 
order and polarization) of spectra can be explained and predicted. The symmetry 
analysis is generic and can easily be applied to other target and/or field configurations. 
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In the past few decades, strong-held physics which studies the interaction between 
strong and ultra-short laser pulses and atoms and molecules has drawn an increasing 
level of interest for the perspectives it offers to probe the organization of matter [U 
[21 S H [5], at the scale of the electronic dynamics in this systems, thus giving birth to 
attosecond science [6j. An efficient means of obtaining such ultra-fast pulses relies on the 
production and manipulation of high-harmonic generation (HHG) emission. For atomic 
targets, our understanding of HHG relies on the conversion of a large number of laser 
photons, previously absorbed by an ionized electron, into a single harmonic photon upon 
recollision [Z11H1I9]. Beyond atoms, similar HHG have been reported and studied for 
molecular targets where one can take advantage of the more complex geometry mm 
compared with the atomic counterpart. 

Numerical investigations of molecular circularly polarized (GP) HHG have been 
reported for specihc geometrical (stretched molecules) and/or laser held (combination 
of elliptic/GP and static laser helds, carrier envelope phase ...) conhgurations [TT] . 
Recently, the interest in GP harmonics have been revived by experimental work using an 
atomic target and a counter rotating bicircular laser [12], opening the way for dichroism 
studies. Using electric dipole selection rules (EDSR), the properties of such bicircular 
atomic spectra have been explained in a photon picture [131HH US]. Here we investigate 
the process of molecular HHG with bicircular laser helds. For molecular systems the 
stakes are high given that EDSR do not apply (because of the absence of rotational 
symmetry) while the more complex geometry ohers perspectives for richer features 
compared to atoms mm- We use discrete time-dependent symmetry (DTDS) [I6] 
analysis and show that molecular GP HHG can be obtained in generic conhgurations of 
symmetric molecules, as shown in hgure [H 

From the theoretical and numerical points of view, most HHG models disregard the 
photon quantum nature of light and compute spectra from the dipole (acceleration) of 
the electronic wavefunction [19] 


Ri 


•HHG 
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d(t) (u) = J^[{^p{t)\a{t)\^l^{t))]{u), 


( 1 ) 


where denotes the Fourier transform and we use the bracket notation. Here {ip {t)) 
is the numerical solution of the time-dependent Schrodinger equation (TDSE) for an 
isolated single active electron (SAE) model taken in the dipole approximation. In the 
length gauge and using atomic units (unless otherwise specihed), the corresponding 
Hamiltonian operator reads 

7i(x,f) =-y+ V(x)-fE(f)-x, (2) 

hoix) 


where Ho is the free-held molecular Hamiltonian, V is the ehective (SAE) potential [T7] 
and E is the laser electric held. 

We present the DTDS analysis for a coplanar bicircular laser held with respective 
frequencies cji and UJ 2 , which has drawn experimental and theoretical interest in the 
past PEI US [IS] but it can easily be generalized to other laser conhgurations. We take 
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Figure 1. Upper part: Equilateral molecular ion (dark green balls) and counter¬ 
rotating bicircular a;/2a; laser field (light gray curves - the arrows indicate the temporal 
evolution of the field). Lower part: Associated HHG spectrum (dark blue curves, left 
hand axis) and polarization (light red curves, right hand axis) for 2 dimensional (solid 
curves) and 3 dimensional (dashed curves) models [13 [T^ . The predicted harmonic 
orders and polarization (from the symmetry analysis) are labeled with the light green 
stripes and upper/lower triangles respectively. As illustrated in the upper part, the 
molecule and laser have compatible symmetries such that, each 1/3 period of the 
laser, the emission signal is repeated with 27r/3 rotation. The combination of all those 
emissions leads to specific circularly polarized harmonics, as shown in the lower part. 


the counter-rotating case as a reference for the laser such that, in the (x, y) polarization 
plane, the bicircular electric held reads 


E{t)=Eof{t) 


coscuit -f cosa;2t 
sin uit — sin u; 2 t 


(3) 


where Eq and / are the (common) peak held amplitude and envelope respectiveljljl. In 
most cases, the beating between the two color components induces a periodic vanishing 
of the laser held (see the light gray curve in the upper part of hgure H]) which is naturally 
associated with recollision and HHG. This becomes obvious in a rotating frame at 
frequency Au = {oji — uj 2 ) /2, which maps Hamiltonian (|2]) to [20] 


Pi (x, t) = T-Lq (x, t) — AuCz + 2i?o/ {t) cos utx, (4) 

where tildes stand for rotating frame coordinates, = —i [xdy — yd^ is the angular 
momentum operator and u = (cji + 0 : 2 ) /2. From the lower part of hgure [3 we notice 
that only some harmonics, with a well dehned polarization, are present in the spectrum 
as we explain in what follows using DTDS analysis. 

The role of dynamical symmetries in HHG computations has been identihed 
previously considering eigenstate decomposition of the wavefunction using free-held [2T] 
or Floquet [TG] |22] state. Here we present a fully non-perturbative analysis of HHG and 


f Absolute phases for each color component have been canceled with an appropriate time shift and 
frame rotation. 
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begin with the laser component of Hamiltonian ([2]). For the theoretical analysis, and 
without loss of generality, we assume 0 < |a;i| < \uj 2 \ in equation ([3]). Finally, in order 
to impose some periodicity on the system, we assume the two color components locked, 
i.e., UJ 2 I 0 J 1 = ^ 2/^1 with ki and k 2 nonzero integers (^2 > 0 thus fci + ^2 > 0) and 
co-prime (irreducible fractional form), which is not a strong hypothesis since rational 
numbers are dense in real numbers M. We illustrate the laser held in counter-rotating 
uj/2uj {k 2 = 2ki) with a light gray curve in the upper part of hgure [H Putting aside 
the envelope (/ = 1 although numerical simulations of the TDSE fully account for it), 
basic trigonometric algebra shows that the electric held (jS]) is 2'Kki/uji-time periodic 
and exhibits a. ki + fc 2 -fold rotational symmetry (for ki + k 2 ^ 0) 

E + ) = ri(n0)E(f), g = 27r , , (5) 

V 1^1 / ki + k 2 


where fl {a) is the rotation matrix around the z axis with angle a, for all integer n. 
For example, in the counter-rotating ujj2u conhguration, this leads to an overall 3-fold 
symmetry (again, see the light gray curve in the upper part of hgure [1]). 

Moving to the molecular part of Hamiltonian (|2]) we assume that the target also 
presents some degree of symmetry which is compatible with the held own symmetry. 
More specihcally, we assume that the potential SAE (ehective) potential presents a k- 
fold symmetry T-Lq (x) = 'Hq (fi {—27r/k)x.). Then, the overall system “molecule-|-laser 
held” presents a fc-fold symmetry, where k as the greatest common divisor between k 
and ki + k 2 

( X, f-I-n-^ J = (O (—n6*) X, f), 9 = 27t^, (6) 

V '^ 1 / k 


as deduced from ([5]), for all integer n. As can be seen in the upper part of hgure [H the 
combination of the equilateral Hj]" molecular ion HU and counter-rotating u:/2u: laser 
held is an example of such a compatible symmetry. Experimentally and numerically, 
HHG corresponds to a range of intensities where ionization is moderate as HHG vanishes 
for depletion of the ground state. Therefore it is reasonable to assume that the 
wavefunction repeats itself over time, following a similar symmetry to equation ([6]) 
- up to a phase that disappears in the dipole (acceleration) computation - and in 
turn, so does the dipole (acceleration). We have checked numerically [I^ the limiting 
factor of ionization by observing that the properties of HHG spectra described in what 
follows degrade when entering a (strong) ionization regime, due to the loss of periodicity 
associated with ground state depletion. 

We combine the symmetry condition (El) with the dehnition of the Fourier 
transform ([1]) over one period of the laser held 
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t + kdL 

uil 

t + k-^' 

UJl 


^-iut 


dt, 


where d^ and dy are the dipole components in the polarization plane and we have used 
simple time translation in each component of the sum in the hnal equality. Then, using 
the symmetry condition ([6]) for the dipole (acceleration) we factorize the spectrum 


R-hhg — 


1 

k 



k=0 


COS kO — sin kO 
sin kO cos kd 



V dy {t) J 

ROe^+ROsy 


^-iut 


dt^ 


where and are the units vectors in the x and ^/-directions respectively and we 
notice that the (constant) factors and Ry have been factorized out of the previous 
sum. Then, using the exponential dehnition of the cosine and sine function, and after 
factorization, the HHG spectrum summarizes to 


R-hhg — 


1 
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2 ^ 
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and the geometric sums can be computed analytically leading to 


(7) 
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k otherwise 


for 


»(l±^ U27rZ 


( 8 ) 


where Z is the set of integers. 

We begin the DTDS analysis with the generic symmetry case k > 2, e.g., for the 
Hg molecular ion and counter-rotating u}/2uj laser conhguration as illustrated in the 
upper panel of hgure [U Since the two frequencies ui and U 2 are mode locked, one 
can dehne a common fundamental frequency u = Ui/ki = 0 : 2 /k 2 and from the overall 
“molecule-|-laser held” symmetry order we dehne the integer h such that ki + k 2 = nk. 
Already, from the trigonometric sums ([8]) we see that all the harmonics 'Lu\k'Luj ± uji 
vanish from the spectrum. For the other ones, with fcZcJ -|- uji the spectrum ([7]) is such 
that Rx/Ry = i which is characteristic of left (counterclockwise) CP emission wile with 
kliUJ — Ui we have Rx/Ry = —i corresponding to right (clockwise) Chl^; in other words, 
the HHG spectrum is circularly polarized with alternating helicities. Going back to our 
molecular model {k = 3 and h = 1) we predict that (3n -|- l)a; and (3n — 1 ) 0 ; are left 
and right CP respectively while 3nu harmonics are absent from the spectrum. Looking 
at the lower part of hgure [H, we have a numerical conhrmation of these predictions for 
two (continuous lines) and three (dash) dimensional simulations. 

The symmetry of the molecular system and laser held do not always match and may 
cancel each other, leading to fc = 1. This is for instance the case for the H^ molecular 
ion with counter-rotating u/3u laser held (fc = 3 and ki + k 2 = 4: - see upper part 


§ We assume oji > 0, otherwise the respective helicities for the khuj ± uji have to the inverted. 










Circularly Polarized Molecular HHG 


6 



Figure 2. Same as figure [T] for the H;^ molecular ion with counter-rotating 
bicircular w/Sw laser field (upper part), with ijj/2uj (middle) and Hj with w/Sw 
(lower) [HIHI]. 


of figure |2]) or H 2 with u/2u {k = 2 and /ci -|- ^2 = 3 - middle panel of the hgure). 
With fc = 1 we see that the sums in equation ([H]) are reduced to a single term and 
destructive/constructive interference effects are lost. As a consequence, all harmonics 
are made available to the HHG spectrum with arbitrary polarization. Here again, this 
is conhrmed numerically in the upper right part of hgure [2] where all harmonics are 
present in the spectrum. 

Finally, among all symmetry orders, the conhguration k = 2 offers an intermediate 
level of prediction. This minimum symmetry order allows constructive and destructive 
interferences in the sums ([8]) leading to only odd harmonics of the common fundamental 
frequency (2Z + 1) cJ to be allowed in the spectrum. However, compared to the generic 
symmetry condition discussed previously, the two subsets 2Zu + Ui and 2ZcJ — Ui 
(respectively left and right CP) lead to the same harmonics in the overall spectrum 
which, in general, prevents drawing any conclusions on the polarization of those 
harmonics. With the Hj molecular ion model, the minimum symmetry conhguration 
is obtained, e.g., with counter-rotating u/?>uj CP helds, as shown in the lower left part 
of hgure [21 Then looking at the associated HHG spectra (lower right part), we have 
the conhrmation that only odd harmonics appear in the spectrum. Note that time- 
dependent or rotational symmetries are not the only discrete symmetries one can use 
to draw predictions on the HHG spectrum. For instance, with a linearly polarized 
laser helds and linear molecules aligned parallel or orthogonal to the laser polarization, 
the system presents a symmetry by rehexion about the polarization axis. In turn, this 
leads to vanishing dipole moment in that specihc direction resulting in linearly polarized 
harmonics (sharing the same polarization direction as the driving laser). 

In conclusion, we have shown that molecular systems irradiated with a bicircular 
laser held oher a very robust framework for the generation of selected circularly polarized 
(CP) harmonics, provided their symmetry order is compatible with that of the laser. 
Using the non-perturbative discrete time-dependent symmetry (DTDS) analysis, we have 
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shown how these harmonics, and their polarization, can easily be predicted. Beyond 
the scope of the present manuscript, the DTDS is relatively generic and can easily be 
extended to other target and/or laser conhgurations. For example, atomic systems can 
be seen as a molecule with vanishing internuclear distance, i.e., the united atom limit. 
In this situation, the total symmetry of the “atom+laser held” always matches that of 
the laser alone because of the rotational symmetry of the atomic target. Then, applying 
a symmetry analysis almost identical to the one described previously we recover the 
predictions of the (perturbative) photon picture where only harmonics (n ± 1) cui + nu 2 
are allowed in the HHG spectrum [131 HH], with opposite CP helicities, for ki + k 2 > 2. 
Note though that in the case of a single CP held (fci + /c 2 = 0) the photon picture 
forbids the generation of harmonics altogether while the symmetry analysis imposes no 
restrictions leading to possibly all harmonics, as shown in [23]. 
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